Stimulated Raman backscattering of laser radiation in deep plasma channels 



^1" 
O 

o 

(N 
>V 



(N 



r/3 
C/3 

o 

C/3 . 



> 



m ■ 

o ■ 
^t- : 
o . 

1/3 ' 

o : 

t/3 1 
Or 



X 



Serguei Kalmyko\Q and Gennady Shvets 
TTie University of Texas at Austin, The Department of Physics, 
Institute for Fusion Studies, One University Station C1500, Austin, Texas 78712 

(Dated: February 2, 2008) 

Stimulated Raman backscattering (RBS) of intense laser radiation confined by a single-mode 
plasma channel with a radial variation of plasma frequency greater than a homogeneous-plasma 
RBS bandwidth is characterized by a strong transverse localization of resonantly-driven electron 
plasma waves (EPW). The EPW localization reduces the peak growth rate of RBS and increases the 
amplification bandwidth. The continuum of non-bound modes of backscattered radiation shrinks the 
transverse field profile in a channel and increases the RBS growth rate. Solution of the initial- value 
problem shows that an electromagnetic pulse amplified by the RBS in the single-mode deep plasma 
channel has a group velocity higher than in the case of homogeneous-plasma Raman amplification. 
Implications to the design of an RBS pulse compressor in a plasma channel are discussed. 

PACS numbers: 52.35 Mw, 52.40 Nk, 52.50 Jm 



I. INTRODUCTION 

Stimulated Raman backscattering (RBS) of laser radi- 
ation in plasmas 0] is a parametric process in which a 
laser beam (pump wave) is backscattered off the electron 
plasma density fluctuations. These density perturbations 
are driven and amplified by the ponderomotive beatwave 
of pump and scattered electromagnetic (EM) waves. Un- 
der certain phase matching conditions a positive feedback 
loop develops that results in the onset of a temporal or 
spatio-temporal instability Q|. The RBS in transversely 
homogeneousplasmas has been extensively studied since 
early 1970s [J, Q], when it first came to the fore in the 
context of fast electron generation and target pre-heat 
in laser confinement fusion. The basic treatment of the 
RBS in homogeneous plasmas is now a classic that can 
be found in a number of textbooks Q . 

Although RBS of long laser beams has been stud- 
ied for at least three decades, the short-pulse regimes 
of this instability have only recently attracted atten- 
tion due to the advances in generation and amplifica- 
tion of sub-picosecond multi-terawatt laser pulses 0. 
The RBS of such pulses in rarefied homogeneous plas- 
mas (cjp -C wo, where luq is a fundamental frequency of 
laser, ui p = \J Ai:e 2 no/m e is an electron Langmiur plasma 
frequency, — |e| is an electron charge, and no is an elec- 
tron plasma, density) was explored in detail in both ex- 
M U (1 S 113 and theory 13 13 13 
The recent upsurge of interest in the RBS 
has been specifically caused by theoretical discovery of 
the possibility to amplify and compress ultra-short laser 
pulses in a plasma by backscatte ring a long low-intensity 
counter-propagating laser beam |l8t lTii |20| . 

Certain applications of short-pulse lasers, such as novel 
X-ray source development |2ll|. generation of high har- 
monics of laser radiation [23] and laser particle accelera- 
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tion |23J,|2J] in tenuous plasmas, benefit from a long in- 
teraction distance. In a homogeneous plasma the region 
of high intensity interaction is confined to approximately 
one Rayleigh diffraction length Zr — k^r^jl, where ko is 
a wavenumber, and ro is a radius of a focal spot of a laser 
beam. Propagation over longer distances requires some 
form of laser guiding. Guiding by a plasma channel is the 
most promising experimental approach [25l 
for high laser intensities. Excitation of relativistic plasma 
waves in channels was analyzed [3 13 13 13] for parti- 
cle acceleration. 

The guided laser pulses arc not immune to paramet- 
ric instabilities, such as forward and near-forward stim- 
ulated Raman scatterin g (S RS) in parabolic 0^ 13 13 
l3rj I37l.l38l.l39t -tapered |4q . leaky |4l| , and single-mode 
flat channels [ij, or large-angle SRS in plasma- filled 
capillaries ji^. Commonly, the energy losses of a laser 
pulse and excessive plasma heating due to the large-angle 
SRS (including RBS) are undesirable for some applica- 
tions |44| , and uncovering new physical mechanisms that 
enable the RBS suppression is up to date. For some ap- 
plications, however, the RBS can be useful. E.g., the 
laser pulse leading front depletion by the RBS may seed 
either forward SRS or resonant modulational instabil- 
ity [3 Novel schemes of short-pulse amplification in 
a plasma[3mi are based on the backscattering of a long 
moderately intense pump laser beam into a short counter- 
propagating signal pulse: the energy of the pump is ab- 
sorbed by the signal as it is amplified and compressed. 
The Raman compression could be a viable path to ob- 
taining high-power single-cycle pulses. For transversely 
homogeneous plasmas, one of the challenges of Raman 
amplification is to ensure the uniformity of plasma along 
the interaction axis so that the signal be downshifted 
from the pump by almost exactly the plasma frequency 
Lu p . In the present paper we discover novel features of the 
RBS in plasma channels, which are favorable for realiza- 
tion of Raman amplification (to anticipate, the RBS pro- 
cess in a deep plasma channel can be broadband enough 
to make exact resonant detuning of pump and signal un- 
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necessary). 

To enable analytic progress and to facilitate qualitative 
understanding, we focus on plasma channels that support 
a single confined hf EM mode, referred to as a fundamen- 
tal mode of channel (FMC) . The laser pulse confined in a 
channel is assumed to be the FMC. The RBS in a single- 
mode channel proceeds differently than in a homogeneous 
plasma and can be characterized by the following novel 
features: (a) reduction of the RBS peak temporal growth 
rate; (b) broadening of the RBS amplification bandwidth; 
(c) modification of the transverse profile of the scattered 
mode from that of FMC. Depending on the parameters 
of channel (such as on-axis plasma density and density 
depression) and pump laser (such as frequency and in- 
tensity), those modifications can be cither more or less 
prominent. 

It was suggested earlier (4^ in the context of Raman 
forward scattering (RFS) that plasma frequency varia- 
tion across the channel may significantly reduce the peak 
growth rate. The effect is resulted from a strong lo- 
calization of a scattering electron plasma wave (EPW) 
near the channel axis. However, due to the compli- 
cated hybrid EM nature of a relativistic EPW in a chan- 
nel [3(], l3l [ I32L I33I I - only approximate results were ob- 
tained in [43. Specifically, it was assumed that the scat- 
tered radiation was in the fundamental channel mode. 
For the RBS, the EM component of the short-wavelength 
EPW may be neglected, enabling us to account for the 
transverse profile modification of the backscattered radi- 
ation. This is accomplished by expanding the transverse 
profile of the scattered field into the channel eigenmodes, 
i.e., the FMC plus continuum of the non-bound hf EM 
modes. These continuum modes of channel (CMC) do 
not exponentially decay outside of the channel (as the 
FMC does) and exhibit a cosine-like behavior at infinity, 
which makes them similar to the transverse Fourier har- 
monics of scattered radiation in a homogeneous plasma. 
The problem is complicated by the fact that the FMC 
and CMC are not independent. Radial shear of the 
plasma density couples them to each other not only cre- 
ating the field structure different from that of FMC but 
also affecting spectral features of the instability. One of 
our goals is to evaluate the continuum-mode contribu- 
tion to the RBS growth rate. We put emphasis on the 
regime of strong plasma wave localization (SPWL) . Spec- 
tral features of this regime are markedly different from 
those of the homogeneous-plasma RBS. The SPWL oc- 
curs when the parameter r\ = (Acj p /wpi) 2 /(27h m/^pi) 
is large compared to unity {here and elsewhere, 7hom = 
y / \ao\ 2 ujo<jJpi/4: is the maximum increment of the weakly 
coupled RBS in a homogeneous plasma Q, Q{, uj p \ = 
lu p (x = 0) is a plasma frequency on axis, Aw^ is the chan- 
nel depth, and ao = eEo / (m e woc) is a normalized ampli- 
tude of electric field of the pump wave}, which physically 
means that the plasma frequency depression in a channel 
exceeds the RBS bandwidth in a homogeneous plasma. 
In the SPWL regime, the maximum temporal increment 
is shown to scale as 7hom/\/?/- Hence, plasma wave lo- 



calization suppresses the instability. On the other hand, 
the transverse shear of the plasma frequency results in 
a broader instability bandwidth. Spectral maximum of 
the backscattered light is found to be red-shifted from the 
pump frequency by more than lo p \ , and the spectrum it- 
self extends on the red side far beyond the frequency 
bandwidth of the homogeneous-plasma RBS. The band- 
width increase is due to the backscattering off the channel 
regions with higher local plasma frequencies and, there- 
fore, higher red-shifts of scattered light. In the SPWL 
regime, the CMC contribution transversely shrinks the 
fastest-growing mode of scattered radiation and concen- 
trates the scattered field near the channel axis. This 
effect is followed by enhancement of the RBS (in the nu- 
merical examples presented in the paper, the continuum 
modes add from 25% to 50% to the peak value of tem- 
poral increment). 

The CMC contribution to the RBS process may be 
neglected if the pump amplitude is sufficiently small 
(or the channel is deep) to satisfy the inequality Uo <C 
(Auip/ujpi) 2 1 (ujo/lu p i) j ^ 4 . This regime is referred to as 
single-mode, because the scattered field is fairly well rep- 
resented solely by FMC. That only three unstable waves 
(two bound EM modes and a localized EPW) participate 
the single-mode RBS makes it the channel analog of the 
three- wave (i.e., weakly coupled) RBS in a homogeneous 
plasma despite the vast difference in spectral features. In 
the single-mode regime an analytic solution of the initial- 
value problem has been found, which describes the evo- 
lution of the single-mode backscattered EM signal in the 
field of a single-mode pump. It is shown that the max- 
imum of the wave packet moves with the velocity 2c/3, 
which is higher than a group velocity of radiation c/2 of 
the weakly coupled RBS in a homogeneous plasma. High 
group velocity of the amplified pulse and its broad band- 
width produced by the transverse shear of plasma density 
profitably distinguish the single-mode Raman amplifica- 
tion in a plasma channel from its homogeneous-plasma 
counterpart [19j. 

The paper is organized as follows. In Sec.^the single- 
mode planar plasma channel is introduced, and basic 
equations governing the nonlinear plasma response to 
the combined pump and scattered radiation are derived. 
These equations are solved in Sec. lIIII bv the mode expan- 
sion, and the generalized dispersion relation is derived. 
This dispersion relation allows for the coupling between 
the FMC and CMC of scattered radiation. In Sec. |IV| 
the dispersion equation is solved in the most interesting 
and novel regime of SPWL. The criterium which allows 
to neglect the CMC contribution to the RBS process is 
proposed. The spectral bandwidth and peak temporal 
growth rate of the single-mode regime are derived. In 
Sec^the initial value problem describing the linear evo- 
lution of the EM signal is formulated and solved, and the 
group velocity of the signal is evaluated. Summary of the 
results is given in Sec. IVII The Appendix reviews some 
properties of the associated Legendre functions with an 
imaginary order which describe the continuous spectrum 
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of the EM channel eigenmodes. 



II. BASIC EQUATIONS 

To begin with, we define the unperturbed state of a 
plasma and laser radiation, that is, the background on 
which the instability grows. In order to make analytic 
progress, a planar plasma channel is chosen with the den- 
sity profile given by 



w 2 (x) = oj 2 p2 - (A^/2) cosh- 2 (x/a). 



(1) 



Plasma frequency in the channel varies between lo p i — 
(uj 2 2 ~ in the center and uj p2 at infinity. Rel- 

ative plasma density depression is n e (x = 0)/n e (|x| = 
oo) = (1 + Ao) 2 /2) _1 , where Aui p = Auj p /uj p i. Through 
the rest of the paper, electron density at the channel 
axis is held fixed, and the channel depth is varied. The 
normalized hf electric field of pump radiation with an 
arbitrary polarization is 



a.o(x, z, t) = Re [ao(:r)< 



ikoz— icuotl 



(2) 



The laser electric field in the channel is the solution of 
the eigenvalue problem 



A)a = [~d 2 /dx 2 + k 2 (x)] a = A a 



(3) 



with the boundary condition ao(|a;| — > oo) — * [here, 
A = (t^o/c) 2 - fc 2 , k p (x) ee u p (x)/c, and a = |a |]. 
In Eq. @, we assume k p \o ~ O(l), k p \ = ui p i/c, and 
ao < 1 and neglect a relativistic correction to the mass 
of electron oscillating in the pump field. Thus, relativis- 
tic self-focusing of a laser beam ^(| is excluded. Rel- 
ativistic self-guiding effects will be addressed in future 
publications. We require that the eigenvalue problem (J3J 
has the unique solution decaying at \x\ — > oo, which gives 
the transverse profile of the FMC 



ao(x) — uq cosh (%/cr) = 



u i>o 



(4) 



and a relation between the plasma channel depth and 
width 



Aw p (t/c = 2. 



(5) 



Also, the eigenvalue equation gives the dispersion relation 
for the pump field 



where fc 



;j2 



\ 2/2 i 2 i 2 

A = uj /c - k = k p2 - cr 
: u p2 /c. 



(6) 



The perturbed hf electric field in a plasma is 

a(aj, z, t) = a + Re [a s (x, z, t ) e ~ iujA+lk ^] , (7) 

where a s = eE s /(m e woc) (|a s | <§; «q) is a complex am- 
plitude of the normalized electric field of backscattered 
radiation, uj s — ojq — uj p i, k s = —ko + k p \. In the case of 
rarefied plasma (u p <C ojq), which is considered here, the 



RBS is a resonant process in which only the Stokes 
component of scattered radiation is involved [see Eq. Q]. 
The amplitudes Oo( s ) are slowly varying in time and space 
on the scales Uq 1 and k^ 1 , respectively. We shall con- 
sider the weakly coupled RBS pi 12; whose temporal in- 
crement is smaller than the electron Langmuir frequency 
(which is valid at uq < ^Juj p /ujq). Hence, the envelope 
of scattered radiation is slowly varying in time and in 
the direction of propagation z: \da s /dt\ <C w p i{2) |a s | , 
\da s /dz\ < k pl( 2)\a s \. 

We neglect the ion density perturbation produced by 
the laser and scattered radiation. We consider ions to be 
a fixed neutralizing positive background in the form of 
a channel with a density profile given by Eq. (I}. This 
assumption is adequate for laser pulses shorter than an 
ion plasma period (rt <§C 27rj~ 1 ). In the opposite limit 
of long pulses, various parametric instabilities have been 
analyzed previously 35] . Ponderomotive force due to the 
interference of incident and scattered radiation excites 
perturbations of electron density, Sh e = n e — n (x), 



6n e (x, z, t) = Re [Sn e (x, z, t ) e ~ l ^ t+lk ' z ] 



(8) 



where k e — fco — k s . In a rarefied plasma, the amplitude 
5n e of the scattering EPW varies in space slowly on the 
scale kQ 1 . Moreover, for the regime of weak coupling, 
this amplitude is slowly varying on plasma temporal and 
spatial periods: \d8n e /dt\ <§; u p i( 2 )\5n e \, \dSn e /dz\ <C 
k p i(2)\5n e \. 

The amplitudes of scattered EM wave and scattering 
EPW obey the coupled-modes equations, which follow 
from the equations of nonrelativistic hydrodynamics of 
electron fluid in the hf EM field Q and Maxwell's equa- 
tions for the scattered radiation, 



d d c 2 

-7T- c— + i— — (£ -A s ) 

at oz 2uj s 



d_ 

at 



-ioj pl {Aoj p /2) 2 (l-Vo) 



a s = UpiGiipov, (9a) 
v = ajpiG 2 4>oas, (9b) 



where v = [uj p (x)/aj p i] 2 [Sn* e /n (x)], Gi = u Q /(4iu> s ), 
G 2 (x) = iw 2 w 2 u with u>o, s = Wo,s/upi, & p = w p (x)/w p i, 
and A s = (lu s /c) 2 — k 2 . The difference between A s and Ao 
is eliminated from all the further equations by the sub- 
stitution a s ,v oc exp[icz(Ao — A s )/(2a> s )]. The set (J^J is 
valid under assumption that the short-wavelength scat- 
tering plasma wave (fc w 2fco 3> fc p i(2)) is predominantly 
electrostatic, with the EM component of the plasma wake 
neglected [42j . We scale time and space to w" 1 and k~i , 
respectively, and introduce the dimensionless variables 
i — LO p it, {z, x} = k p i {z,x} (all the overbarred quanti- 
ties which appear in the equations below are normalized 
in this way). Eqs. 10 recast now as 



9 d i - s 

M~dz + Ws {a ' a} 

|_-z(A^/2) 2 (l-^ 2 ) 



a s = Giipov, (10a) 
v = G 2 tp a s . (10b) 
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In Sec. IIIII we derive from Eqs. (|10fl the dispersion rela- 
tion of RBS in the channel QJ. 



III. GENERAL DISPERSION RELATION 



with s = fi = 1. The degree s = 1 given, the spectrum 
of Eq. 1)12(1 consists of one discrete level with /i — 1 and 



a continuum of modes with fi 



q , q real. Hence, the 



full set of solutions of Eq. is composed of the single 
FMC Pl(y) = tp and the set of CMC 



Eqs. Q1UJ1 can be solved by using the Fourier-Laplace 
transform of the envelopes, a s , v oc exp(i(Di — ikz). It is 
convenient to introduce a new variable y = tanh^/cr), 

-,2 



y 



1 + C 2 y 2 , where C 2 



and express ipo 
A<D 2 /2 = 2 /a 2 . Plasma occupies the area — 1 < y < 1. 
Expressing v through a s from Eq. (|10b|l , inserting it into 
Eq. l|lUa|) . and using Eqs. (J3J, ©, we obtain the equation 
for a s (y, k,u>), 



r -2 

l-\\a. a = a 



y 2 u 



(ii) 
(i + 

The 



where p 2 = 2uj s (uj + k), G — u 2 ,u 2 /2, Q(co,y) 
C 2 y 2 )/(1 - B 2 y 2 \ B 2 = A^ 2 /(4cS) = (a 2 ^ 
function p 2 {k) contains the information about propaga- 
tion of scattered radiation. The operator in the LHS of 
Eq. (|1 If) comes from the equation for the associated Leg- 
endre functions Ha 



P? iq (y) = 



tanh(a;/<r) =F iq 
(l T iq)T(lT*q)' 

i y t n ( i 



r(i Tiq) i Tiq \i-y 



D ±iqxfo 



y 



±iq/2 



(13) 



At | a; | — + oo, the CMC ((13(1 reveal the cosine-like behavior 
[see Eqs. (|A2() ]. Unlike the pump field given in the form 
of the FMC, the scattered radiation in a channel is not 
necessarily represented by FMC only. The CMC 1(13(1 
describe the discrepancy between the true transverse field 
profile and the FMC. Therefore, we express the general 
solution of Eq. 111(1 as a mode expansion 



a s (y) = Pi(y) 



Pl 9 (y)a s (q)dq. 



(14) 



C^(y)={ _(l-y 2 )_ + s (.s + l 



dy 



dy 



>Pj%) = 
(12) 



Knowing that C h iPl = 0, £ M Pf(y) 
y 2 )~ l P\ q {y), we arrive at the equation 



(1+9 2 )(1 



+ OG 



Pl q (y)- 



ja s (q) dq = -a 2 



G 



+ 00 



i -y 2 



+ -g(Q,y)) \Pt(y)+ / P?(y)a s (q)dq 



(15) 



We multiply Eq. 1(15(1 by P\{y) and integrate it over y. where 
Having in mind that J^jP^iy)^ - y 2 )~ 1/2 dy = [for- 
mula (7.132.1) of Ref. 1471. we arrive at the dispersion , .. , , 

F(w,iq)= / Q{oj,y){l-y 2 ) 1/2 Pl q {y)dy, (17) 



equation 



G G f 

p 2 H —Q(^) — — / F(&,q)a>s{q) dq, (16) and the plasma response function reads 

2uj 2lv J 



l 

(l-y 2 )g(u,y)dy=-^(l- — + ^ 



(B 2 + C 2 ){1 - B 2 ) , fl + B 
In ' 



B 5 



1 - B 



(18) 



It has been shown previously 01 that the function l(18|) the limit of a broad shallow channel, that is k e a ^> 1. In 
describes a purely electrostatic response of a plasma in the case of RBS k e ss 2/cq. Thus, any channel with a > 
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k^ 1 is wide, and the electrostatic description is always 
applicable to the plasma response in the RBS process. 

In order to find a closed-form dispersion equation we 
have to determine explicitly the CMC amplitudes a s (g) in 

Eq. JH3). We multiply Eq. JTHJ by Pp 9 ' (y) and integrate 
it over y. Using the orthogonality condition lAlfl and the 
identity T (1 + iq)T (1 — iq) = irq/ sinh(7rg) [4j], we obtain 
the expression 

2 ( G \ q/ sinh(7T<j) 



2Q J 1 + (ap) 2 + q 2 

+00 

F(ui,-iq)+ / Fi(u) : q,q')a s (q)dq 



(19) 



where the kernel iq = f\ G{&, y)P\ q (y)Pi iq {y)dy de- 
scribes the coupling between the modes of continuum 
with different indices q. Further, we take into account 
only the coupling between the FMC and CMC, and drop 
the last term in brackets in the RHS of Eq. (fT^|l . The 
resulting dispersion relation reads 



<fr(u;,p) 



F(u),iq)F(u,-iq) qdq 
1 + {op) 2 + q 2 sinh(7rq) 



(20a) 
(20b) 



Eq. (|20a|l describes the spatio-temporal evolution of ini- 
tial perturbations of field and electron density 0| ■ Real 
and imaginary parts of the complex frequency Q as a 
function of a real wavenumber shift k determine the tem- 
poral evolution of spatial Fourier harmonics of the signal, 
and can be used to solve the initial-value problem. In 
the next Section, the complex solution ui{k) of the dis- 
persion equation H20al) (with k real) will be found in the 
limit of strong plasma wave localization (SPWL). The 
SPWL is achieved when the channel depth Aid 2 is much 



larger than the instability bandwidth, i.e. Acu 2 ^> \u}\, 
or \B\ 2 > 1. The unstable EPW will be found confined 
in the near-axis area with a transverse extent of about 
|-B| _1 . For this regime, the contribution of CMC to the 
growth rate can be important. 



IV. SOLUTION OF DISPERSION RELATION IN 
THE SPWL REGIME 

A. Temporal increment of instability 

The parameter area for the SPWL regime is prescribed 
by the inequality \B\ 2 > 1, or \u>\ < (Aui p /2) 2 . Taking 
the estimate \lo\ ~ 7hom = (uo/2)tJujq, we establish the 
limitation on the laser amplitude from above 

u < U = Aw£/(2V5o). (21) 
Eq. H21fl provides the same scaling in uiq as the condi- 
tion of weak coupling u$ < for the RBS in ho- 
mogeneous plasmas 2J. In the SPWL regime, the unsta- 
ble electron density perturbations are localized stronger 
than the EM waves. The size of the area of localization 
is Sx - 2a/\B\. 

To evaluate the RHS of the dispersion equation l|20a|) 
in the limit \B\ 2 ^S> 1, we note that the integral l|17fl that 
determines the kernel of Eq. (|20b|l can be represented as 
F(Q, iq) = j\ T(y, iq)(l - B 2 y 2 )- 1 dy, where T{y, iq) = 

{l + C 2 y 2 )Pl{y)Pl q {y) is independent of B. At \B\ 2 > 1, 
only y values from the close vicinity of channel axis, that 
is, \y\ < \B\~ X , will contribute to the integral. Therefore, 
effective coupling between the FMC and CMC occurs 
near the channel axis in the area with a transverse extent 
of order In this case, the RHS of Eq. I|20all can 

be expanded in powers of B . To evaluate the lowest- 
order term of the expansion, the on-axis value of T(y, iq) 
is taken, i.e., T(0,iq) = —iq/[(l — iq)T(l — iq)]. Then, 
the integral (|20b|) becomes 



$(|B| 2 >1) 




dy 



+ 00 



1-B 2 y 2 



q 2 dq 



[l + {ap) 2 + q 2 }{l + q 2 ) 



l+(ffp)-l 
(oBp) 2 



hr 



1+B 

l-B 



l + (apf 



{oBp) 2 



(22) 



Evaluating the plasma response function in the limit of 
large \B\ as Q{uj) ~ iir/B, we reduce Eq. (|20a|l to a 
relatively simple form: 

The RHS of this equation represents the contribution 
from the CMC which cannot be a priori neglected. The 



generic SPWL regime can be subdivided into the single- 
and multi-mode sub-regimes. Below, we find the bound- 
ary between them and their spectral features. 

Under the condition \ap\ 2 - 2ui /\B\ 2 < 1, the CMC 
contribution may be eliminated from Eq. 123|) : 

p 2 + [G/(2uj)](in/B) =0. (24) 

Eq. (|24|l describes the interaction of the pump wave 
(FMC) and the single FMC of scattered light via the 
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strongly localized EPW. Therefore, the single-mode sub- 
regime admits an analogy with the three-wave RBS in 
a homogeneous plasma (these processes, however, have 
quite different spectral properties). Eq. (|24|l yields the 
solution with a maximum imaginary part, 

Q(k = 0) = e™/ 3 V^Wv) Thom, (25) 

where the parameter 77 = A(Dp/(27h om ) > 1 is of the 
order of \B\ 2 . Real part of the solution l|25|) gives the red- 
shift of the spectral maximum Aw = \/n 2 / (2r])(% orn /2). 
Eq. I|24|) also predicts a blue-side limitation of the RBS 
bandwidth at Re<I> c = — ^/n 2 /r]( ; yh om /2). From the red 
side (Re a; > 0), Eq. gives no limitation, and the 
spectrum has a tail far extended in this area. 

The validity condition \B\ 2 > 2d> > 1 for Eq. (23} sep- 
arates the SPWL sub-regimes. It is more restrictive than 
\B\ 2 > 1 and may be written as \u)\ < Aid 2 / (8u> ) . Sub- 
stituting the solution l|25(l into the latter inequality pro- 
vides the parameter area for the single-mode sub-regime: 

u Q < U x = C/ /[V^(2cD ) 3 / 4 ] « U . (26) 

The larger ito, the stronger the instability is driven, and 
the larger is the population of CMC excited. When uo 
falls within the interval 

f/i < u < Uo, (27) 

the CMC contribution is no more negligible, and the 
RBS process becomes essentially multi-mode. The RHS 
of Eq. (|23|l cannot be omitted then, and the dispersion 
equation i|23|) is solved numerically. 

In Fig. ^ we present an example of the dispersion 
curves obtained for the RBS in both single- and multi- 
mode SPWL regimes. For the fixed values of normalized 
laser frequency <I>o = 10 and amplitude of electric field 
i*o = 0.007 the increment Imw is found numerically ver- 
sus Rea; for two different values of the channel depth: 
(I) A(D p = 2, or n e (x = 0)/n e (\x\ = 00) = 1/3, and (II) 
AQ p = 1/V2, or n e (x = 0)/n e (|a;| = 00) = 0.8. In Fig.HJ 
the solution of full Eq. Il2.'ill for the set of parameters (I) 
is plotted with the curve (1). The curves (2) and (3) are 
obtained via numerical solution of Eq. H23fl for the set of 
parameters (II), with the CMC contribution deducted in 
the case of the curve (2). The reference spectrum of RBS 
in a homogeneous plasma is presented in Fig. ^ by the 
curve (4). 

Parameter sets (I) and (II) correspond to the SPWL 
regime as the condition uq <C Uq is very well satisfied 
for both. For the set (I) the inequality uq < U\ ~ 0.019 
holds, and the single-mode regime is the case. Contribu- 
tion of the CMC is negligibly small: numerical solutions 
of the full [Eq. (|23J)] and single- mode [Eq. (P3}] dispersion 
equations coincide within the thickness of the line (1). 
Such a good coincidence is not the case for the parame- 
ter set (II) , which corresponds to the multi-mode regime 
(uq > U\ « 0.0047). Comparison between the curves (2) 
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FIG. 1: Spectra of RBS in the regime of strong plasma wave 
localization. Normalized increment is plotted against the real 
part of frequency. Laser and plasma parameters are u>o — 10 
and uo ~ 0.007. Plasma frequency depression in a channel 
is (I) Aujp = 2 [curve (1)] and (II) Aui p = l/y/2 [curves (2) 
and (3)]. The line (4) is a reference curve corresponding to 
the homogeneous plasma (Gj p = 1). Lateral variation of ui p 
red-shifts the spectral maxima and creates the tails on the 
red side. 

and (3) demonstrates considerable enhancement of RBS 
due to the CMC contribution: the CMC correction to 
the peak increment amounts to about 25%. The basic 
characteristics of RBS in the multi-mode SPWL regime 
can be summarized as follows: 

1. The peak growth rate is reduced if compared with 
the case of homogeneous plasma. 

2. Contribution from the CMC enhances the scatter- 
ing process. 

3. The RBS bandwidth inside a channel is signifi- 
cantly larger than in a homogeneous plasma. 

4. The frequency spectrum experiences overall red- 
shift from the Stokes frequency loq — uj p \ . 

The last two features are clear advantages of the SPWL 
regime for the Raman amplification of short pulses in 
plasmas [l^. Due to the broadband nature of the pro- 
cess, exact tuning the signal frequency to ujq — tu p i is 
not necessary to get a considerable amplification rate in 
the linear regime. We can also estimate here the RBS 
growth rate modification for the parameters typical of 
a channel-guided laser driven accelerator. The plasma 
channel created in the experiment [2{j was capable of 
single-mode guiding of the laser pulse with a radius 8^tm 
at the level e -2 in intensity, which under ansatz gives 
a = 4.65/im. The electron density at the bottom of the 
channel ri\ = 4 x 10 18 cm -3 gives k p i<j = 1.75 and, ac- 
cording to the matching condition JSJ), provides the effec- 
tive normalized electron density depression Auj 2 w 1.3 
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FIG. 2: Transverse profile of intensity for the fastest growing 
spectral component of scattered radiation (solid line) numer- 
ically evaluated from Eq. 1(14(1 for the parameter set (II) of 
Fig. complex frequency given by the maximum of the curve 
(3). The dashed line is an intensity profile of the bound mode. 
CMC contribution shrinks the transverse profile of scattered 
radiation and increases the field near the channel axis. 



(we have to mention that an actual channel shape [29| 
is a plasma column with a density depression at the 
axis surrounded by the walls which are few laser wave- 
lengths thick and approximately twice the bottom den- 
sity; as the walls density is much below critical, we as- 
sume that the main contribution to the guiding effect is 
made from the near-axis density profile which admits the 
approximation in the form Q with the effective channel 
depth just calculated). The laser wavelength Ao = 0.8/im 
gives loq w 20. To fall within the SPWL regime of 
RBS the guided laser pulse must possess the intensity 
/ < 4 x 10 16 W/cm 2 to be in accordance with the limi- 
tation M < U w 0.14 [sec Eq. l|2T))], For J«3x 10 14 
W/cm 2 (tto « 0.012, 7 hom w 0.028w p i) RBS proceeds 
in the multi-mode SPWL regime with rj w 24. The in- 
crement formally evaluated from Eq. 1(25(1 gives the peak 
increment 7 w 0.57h O m whereas taking account of the 
CMC contribution increases it to 0.87h om - Therefore, in 
the regime considered, 20% reduction of the peak RBS 
increment can be expected in comparison with the case 
of homogeneous plasma, and the contribution from CMC 
to the peak increment amounts to 30% of Thom- 



B. Transverse profile of scattered radiation in 
multi-mode SPWL regime 

Results of the previous Subsection show high sensitiv- 
ity of the RBS to the transverse structure of scattered 
radiation in a channel: taking account of CMC increases 
the value of the growth rate. Here, we find the CMC- 
related correction to the transverse profile of scattered 



radiation. For the parameters of the case (II) [the curve 
(3) in Fig.^ we evaluate the integral 1)14(1 numerically for 
the fastest growing mode, and present the transverse dis- 
tributions of intensity |a s (a;)| 2 in Fig- EI T ne solid curve 
is given by Eq. 1|14|) . The reference intensity profile of the 
FMC is plotted with a dashed line. The coupling between 
the FMC and CMC transversely compresses the scattered 
radiation beam. To evaluate the effective compression 
numerically we define the rms beam size as y (x 2 ) = 

I +00 +00 x 1 / 2 

[f-™x 2 \as{x)\ 2 dx/ \a s (x)\ 2 dx) . For the inten- 
sity profiles with and without CMC shown in Fig. [5] 
the ratio of the rms sizes is J (5fmc+cmc) / (^fmc) ~ 
0.77. So, the numerical example shows 23% trans- 
verse compression of the intensity profile reached in 
the multi-mode regime. Simultaneously for the pa- 
rameters chosen, power of the scattered radiation re- 
mains almost unchanged: difference between the inte- 
grals /+" |a s (x)| 2 dx calculated for the FMC and FMC 
+ CMC is about 3%. Hence, in the multi-mode SPWL 
regime, coupling between FMC and CMC compresses the 
scattered radiation near channel axis where the pump 
field has maximum. As a consequence, increase in the 
peak increment occurs. For the estimates made at the 
end of previous Subsection the compression was about 
30% and intensity enhancement on the axis about 40%. 



V. SPATIO-TEMPORAL EVOLUTION OF 
BACKSCATTERED PULSE IN THE 
SINGLE-MODE SPWL REGIME 

A. Group velocity of scattered radiation 

In this Section we address the linear Raman amplifi- 
cation of EM wave packet in a plasma channel. First, 
we evaluate the group velocity of backscattered light, 
i.e., v g = c(dRecD/dk) at dlmui/dk = 0. By definition, 
thus calculated group velocity determines the speed of 
the pulse peak. On substituting Q = fl + with fl and 
7 real into the dispersion relation 1)24(1 ■ we separate real 
and imaginary parts of the equation, exclude 7, and ob- 
tain the algebraic relation which defines implicitly the 
dispersion function £!(/c): 

k 2 + Akn + m 2 -~ — — = 28 

3ft + 2k y ' 

Differentiating Eq. I(28|) with respect to k, and plugging 
into the resulting equation k = and f2(0) = Reu) from 
Eq. ((25(1 (i.e., the wave number and real part of frequency 
corresponding to the peak growth rate), we find the group 
velocity v g = — 2c/3 (the minus sign means that the am- 
plified pulse moves in the backward direction). 

The absolute value of the group velocity in a chan- 
nel is higher than in a homogeneous plasma, where 
v ghom — — c/2. To get a qualitative interpretation of 
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this effect we return to the basic equations l|10|) . Tak- 
ing the single-mode approximation for the scattered light 
envelope a s (x,z,t) ~ A a (z,t)ipo(x), initial condition 
u(t — 0) = 0, and the plasma response function Q(u>) « 
m/B{u)) = (2tt/ Au! p )y/iQ, we reduce the set (fTTTft to 
the single equation (c% — d s )A s « u [Aa>p/(16a)o)](-/V(i)), 
where the overlap integral between the scattering elec- 
tron density perturbation and transverse intensity pro- 
file of the FMC (N(t)) = ip%(6n*Jn Q ) dx = 



uoV^i{2uja I Aujp) 2 J Q A s (r)/\ft — r dr describes the ef- 
fective plasma response. The transverse shear of 
the plasma density produces an effective decay of the 
plasma response expressed in terms of the convolu- 
tion (N(t)) cx A s (t) * t" 1 / 2 [compare with the non- 
damped case of homogeneous plasma, where (dj— dz)a s = 

7hom/o o-s(r)dT = jl om a s {t) * 1]. Therefore, the tail of 
the amplified signal experiences the growth rate reduc- 
tion according to A s (t) cx (N(t)) cx A s (t) * i -1 / 2 , and, 
consequently, the signal maximum moves closer to the 
signal leading front than in a homogeneous-plasma case. 
This argument qualitatively explains the increase in the 
group velocity of scattered light in a channel. 



B. Evolution of EM wave packet 

We solve here the initial-value problem for the equa- 
tions lfTU|l 

a s (x,z,i=0) = a s0 (z)ipo(x), (29a) 
P(x,z,i= 0) = 0. (29b) 

which specifies the EM wave packet initially matched 
with an unperturbed plasma channel. We naturally take 
account of these initial conditions by introducing a new 
dependent variable a\(x,z, t) — a s (x, z,t)H(t) [where 
H(t) is the Heaviside step- function] . In the character- 
istic variables 9 = —z and r] = —z — i, the set QlUfl reads 



+ 7T^(£o-A ) 



2lD 



_c9_ 

d .AlS 
drj ii 



v+G2ipoal = 0, 



where 5(t) is the Dirac delta-function. We ap- 
ply the Fourier transform with respect to the co- 
moving "spatial" variable "q, assuming 5(9 — rj) = 
(27r) _1 J _^ exp[ifc(?7 — 6)]dk, and then exclude the 
Fourier image v(x, 9, k). Multiplying the resulting equa- 
tion for the Fourier image a\(x,9, k) by ipo(x) and then 
integrating over x we obtain the averaged equation 



1 



<V>o 2 > 



^#; + 2^ 



(ipo (£ "A ) al) 



(30) 



-G • 



V 3 [i+c 2 (i-^ 2 )] 

fc[l-^ 2 (l-^o 2 )] 



■a aO (-0)e 



-ik0 



where B\ = — Aw 2 /(4fc). In a generic case, trans- 
verse profile of the signal a\(x) is a superposition of 
the FMC and CMC. In this Section, we consider the 
single-mode SPWL regime of Raman amplification and, 
for all instants of time, take the wave packet in the 
form of FMC [hence, (£ — Ao)a| = 0]- We substitute 
a\(x, 0, k) = A(9, k)ipo[x) into Eq. (JSUJl and get finally 



(Q(k)\ 



dA 2 , . , 

de + ^[-2ik) A 



a s0 (-9)e 



-ike 



(31) 



Here, the plasma response function [compare to Eq. IT 
reads 



Q(k) 



1 + C 2 y 



1 



V 2 ) dy- 



(32) 



In the SPWL regime, relatively large wave number shifts, 
\k\ > (AcDp/2) 2 (that is, \B k \ 2 > 1), contribute to the 
signal evolution. Hence, the main contribution to the in- 
tegral (|32J) is made by the integrand values in the vicin- 
ity of y = (or \y\ <C |-Bfc| _1 ). Approximating the inte- 
grand as (l — -B 2 y 2 ) 1 we find the approximate response 
function Q(k) ~ — 2irVk/ Alj p which allows to present 
Eq. fpTTf) in the form 



(d/86- 



iTk 



-1/2 



A = a s0 (- 



ikO 



(33) 



where T — TTThom yTtom/ (277). Solution of Eq. I|33(l reads 



A(6,k) = 1 




^-ikie-e^-i^/Vk d9l _ ( 34) 



Inverting the Fourier transform l|34|l and returning to the 
lab- frame variables, we find the longitudinal evolution 
of the signal, A(z,t) = f °°a s o(6>i + z)V(9 1 - t,9 1 )d9 1 , 
where the the Green function of RBS in the single-mode 



SPWL regime, V{d x -t, X ) = (27T)- 1 J_ ^ exp[ik(dx-t)- 

can be found explicitly in terms of a gener- 
alized hypergeometric function 1F3 |49|: 



V(6 1 -t,6 1 ) = 5(t-6 1 ) + iH(t-0 1 ) (F0-i_y 

l F :] [ll,^l l ^9 1 f{t-9 1 )/^j. (35) 

Solution of the initial- value problem expressed explicitly 
in terms of the Green function (|35|l reads 
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a\(x,z,t) = tp (x) < a sQ (t + z) + iT 2 / iF ; 



a s o( 



9f dO x 



(36) 



-3 

x 10 
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pl 

FIG. 3: Longitudinal on-axis profiles of intensity of the am- 
plified signal for the parameter set (I) of Fig. Q t\ = 1800, 
I2 = 2100, ti — 2400. The seed signal is Gaussian, a s o — 
(mo/10) exp(-z 2 /!00). 

In Fig. |3 we present the evolution of initially Gaussian 
signal a s0 — (uo/10) exp (— z 2 /100) for the parameter set 
(I) of Fig. 1. Snapshots of the intensity profile are pre- 
sented for three consequent instants of time. The front 
of the amplified signal moves with a speed of light in the 
negative z direction, whereas its maximum moves with a 
group velocity v g = — 2c/3, as predicted in the previous 
Subsection. The maximum of the amplified signal grows 
exponentially in time with an increment twice as given 
by (|25|) . Therefore, the solution of initial- value prob- 
lem confirms the predictions of the dispersion analysis of 
Sec. HV| 



VI. CONCLUSION 

We have investigated the RBS of laser radiation in 
the regime of strong plasma wave localization (SPWL) 
in a plane plasma channel, which can support only one 
trapped EM mode (the single- mode channel |42|). For 
the SPWL regime, transverse variation of the plasma fre- 
quency exceeds the RBS bandwidth calculated for a ho- 
mogeneous plasma with an on-axis electron density, and 
the scattering plasma wave is localized stronger than a 
driving beatwave of pump and scattered radiation. In 
this case, the transverse profile of the scattered radia- 



tion is a superposition of the fundamental mode of a 
channel (FMC) and a continuum of non-bound modes 
of channel (CMC). Depending on the plasma and laser 
parameters, the RBS with SPWL can proceed either in 
a single- or multi-mode regime. In the single-mode case, 
excitation of the CMC is suppressed almost completely, 
and only the FMC of scattered radiation is involved in 
the process. This allows for a physical analogy between 
the single-mode SPWL regime of RBS in a plasma chan- 
nel and a three-wave RBS in a homogeneous plasma. In 
the multi-mode SPWL regime, the CMC play an essen- 
tial role. The CMC contribution transversely shrinks the 
scattered radiation beam, and increases the growth rate 
of the instability. Spectral features of both sub-regimes 
are qualitatively similar. The temporal growth rate is 
slower, and amplification bandwidth is greater than in 
the case of homogeneous plasma, and the frequency spec- 
trum experiences an overall red-shift. The group velocity 
of scattered radiation, v g w — 2c/3, is increased versus 
its homogeneous-plasma value, v ghom = —c/2. All these 
features are the consequences of the transverse shear of 
the electron plasma density. The broadband nature of 
RBS in the SPWL regime and relatively high group ve- 
locity of scattered radiation are the features profitable for 
the Raman amplification of short pulses in deep plasma 
channels. 
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APPENDIX A: PROPERTIES OF ASSOCIATED 
LEGENDRE FUNCTIONS WITH IMAGINARY 
ORDER 

The associated Legendre functions P^ %q (y) [H(J , where 
q is a positive real number, are the solutions l|13|) of 
Eq. (|12fl with s = 1 and /1 = ±iq. The orthogonality 
condition for P 1 lq (y) reads 

1 

1 fuii/ \ d y sinh(Trq) , 

7T / p i (y) p i (yh — 2 = % -9), Ai) 

-1 

which is evaluated using the substitution y = 
tanh(x/cr), formulas (3.982), (3.987.1), and (3.981.6) H3 
and the Dirac delta- function property 8{q' — q) = 
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(27T)" 1 /+~ e ix (i-i') dx. The functions Pf^x) do not V\ q {x) = [P[ q {x)+P[ q {-x))/2 andPp 9 ^) = [Pi 19 ( x ) + 
vanish at \x/a\ — > oo and have the asymptotic P 1 _II? (— x)]/2: 

. / r \ piqx/a 

*(*^)~w=«y (A2a) 

^(r- 00 ) ~ l^rJf(T^)- (A2b) C0&((Zx/ff) -i2Ji m r(i + ^ + (Z Rer(i + ^)- 

The cosine function can be constructed of P 1 zq (x) 
by summing up the even complex conjugate solutions 
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